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the  anomalous  reflection  phenomenon  is  derived.  This  theory 
permits  the  determination  of  peak  pressure,  shape  and  duration 
of  the  positive  pressure  pulse  and  appears  to  be  applicable 
to  shallow  explosions  in  deep  water  at  depths  greater  than 
3  charge  radii  from  the  free  surface. 
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The  study  reoorted  here  was  carried  out  under  Task  Re2c-65-l 
as  part  of  a  program  aimed  at  investigating  the  anomalous 
surface  reflection  of  shockwaves  from  shallow  explosions  in 
deep  water.  Since  the  material  in  this  report  is  of  immediate 
interest  to  other  groups,  it  is  published  in  its  present  form 
at  this  time.  The  experimental  data  with  which  this  work  can 
be  compared  will  be  reported  later. 

The  authors  wish  to  point  out  the  essentially  non-stationary 
nature  of  the  entire  problem,  and  the  serious  assumptions 
inherent  in  the  stationary  flow  simplification  U3ed  in  this 
paper.  The  senior  author  (H.  G.  S.)  has  suggested  a  possible 
attack  leading  to  a  more  exact  solution,  which  will  be  pre¬ 
sented  later  if  time  permits. 

This  report  is  for  information  only,  and  the  opinions  ex¬ 
pressed  herein  are  those  of  the  authors. 

This  revision  supersedes  NAVORD  Report  2710,  dated  1  January 
1953. 


W.  W.  WJLBOURNE 
Captain,  USN 

ft.  E.  Ablard 
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ON  THE  OEL.IQUE  REFLECTION  OF  UNDERWATER  SHOCKWAVES 
FROM  A  FREE  SURFACE  I 

I  INTRODUCTION 

When  a  shockwave  produced  by  an  underwater  explosion 
comes  in  contact  with  a  free  surface,  its  pressure  must  drop 
to  approximately  atmospheric,  or  more  precisely,  to  the  same 
pressure  as  the  air  shock  produced  at  the  interface.  As  a 
result,  a  tension  wave  is  reflected  from  the  surface,  which 
modifies  the  incident  wave  in  the  water.  At  the  free  surface, 
the  situation  of  oblique  reflection  of  a  shockwave  can  be 
treated  as  a  quasi-stationary  flow  problem,  where  the  pressure 
is  relieved  by  a  flow  around  a  Prandtl-Meyer  corner  [L]*.  As 
was  first  pointed  out  by  Penney,  the  flow  into  the  expan¬ 

sion  wedge  becomes  subsonic  if  the  glancing  angle  of  the  shock- 
wave  is  sufficiently  small.  As  a  result,  the  rarefaction  wave 
overtakes  the  incident  shock  front  and  attenuates  it.  This 
leads  to  a  region  of  peak  pressure  distortion,  whose  boundaries 
have  been  calculated  by  Penney  and  Keil  [I]  for  a  spherical 
shockwave . 

In  the  present  paper,  an  attempt  is  made  to  analyze  the 
phenomenon  within  the  distorted  region.  A  pseudo-stationary 
approach  is  used,  and  as  a  result  it  appears  that  the  pressure 
relief  in  this  region  takes  place  in  two  stages:  Part  of  the 
pressure  is  released  by  means  of  a  wave  which  is  traveling 
faster  than  the  incident  shock;  the  remaining  pressure  is 
then  relieved  through  a  Prandtl-Meyer  corner  flow.  On  the 
basis  of  this  treatment,  it  is  possible  to  construct  pressure¬ 
time  histories  for  shockwaves  in  the  regions  of  anomalous 
reflection. 


*A11  such  letters  refer  to  the  list  of  references  at  the  end 
of  this  report. 
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II  REGULAR  REFLECTION  AND  THE  PRANDTL-MEYER  CORNER  FLOW 

The  interaction  of  a  plane  shockwave  of  infinite  duration 
with  a  free  surface  at  a  sufficiently  large  glancing  angle, 
can  be  represented  by  a  stationary  flow  pattern  if  a  steady 
flow  is  superposed  parallel  to  the  surface.  This  flow  is 
shown  by  the  vector  -  U  sec  ®  in  Figure  1,  where  U  is  the 
propagation  velocity  of  the  shockwave,  and  ®  is  the  angle 
between  the  wave  front  and  the  normal  to  the  surface.  Since 
a  positive  disturbance  of  finite  amplitude  must  travel  with 
supersonic  velocity  (U>cQ,  the  sound  velocity  in  the  undis¬ 
turbed  medium),  this  superposed  flow  is  always  supersonic. 

This  flow  corresponds  to  the  undisturbed  water  before  the 
shock  front,  as  observed  from  the  moving  point  of  intersection 
of  the  shock  front  with  the  free  surface.  The  flow  behind  the 
shock  front  is  modified  by  a  quantity  u,  the  actual  particle 
velocity  behind  the  front.  This  flow,  f1,  can  be  obtained  by 
vector  addition  of  the  original  superposed  flow  and  u,  as 
shown  at  the  bottom  of  Figure  1. 

The  pressure  at  the  free  water  surface  must  remain  atmos¬ 
pheric.  As  long  as  the  flow  f1  is  supersonic  (or  sonic),  this 
condition  can  be  realized  by  a  stationary  expansion  wave, 
centered  at  the  point  of  Intersection  of  the  shock  front  with 
the  surface.  This  is  the  well-known  Prandtl-Meyer  corner 
flow.  Figure  2  illustrates  this  flow  more  clearly;  the  expan¬ 
sion  wedge  consists  of  the  Mach  lines  of  the  corner  flow. 

According  to  the  Prandtl-Meyer  analysis  of  a  2-dimenslon- 
al  irrotatlonal  steady  supersonic  expansive  flow  around  a 
corner  [L],  the  relationship  between  the  angular  deflection  of 
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FIG.  I  REGULAR  REFLECTION  OF  PLANE  SHOCKWAVE 


FIG.  2  PRANDTL-MEYER  CORNER  FLOW 
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the  Mach  lines,  rj,  and  the  pressure,  P,  is  given  by  the 
expression 


d 


<1  +  If )« 

+  2cd^  -  2co  - 


7W* 


(i) 


c  and  V  designate  the  sound  velocity  and  the  specific  volume, 
while  oo  is  the  specific  enthalpy,  defined  by 


oo  =  jvdP 


(2) 


ooi  is  the  specific  enthalpy  of  the  flow  fi  at  the  beginning 
of  the  expansion. 

Substitution  of  numerical  values  from  the  Appendix 
permits  analytic  integration  of  differential  equation  (1)  over 
the  pressure  range  of  interest.  For  a  complete  expansion  to 
zero  pre^ure,  the  integral 


(*P=0  r  ?  fi 

rj  =  j  dr|  =  |^(l-9.046  *  10"^N-1.048  *  10  P) 


P=0 


(N-2.318*  lO^P)1/2 

Jp=p 


(3) 


is  obtained,  in  which  P,  the  pressure  in  psi,  is  evaluated 
between  the  limits  Pi,(the  pressure  of  the  flow  f^,  and  0. 

N  is  a  function  of  the  flow  f1  into  the  corner  and  is  given  by 

f  2 

N  =  -  1  +  6.282  *  10~6P±^  2.318  *  10'^  (4) 
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The  equality  sign  on  the  right  hand  side  refers  to  sonic 
entry  into  the  corner.  In  the  latter  case, 

*  0.276  Pi1/2  (5) 

for  moderate  pressures. 

As  a  first  approximation,  the  change  in  angle  rj  corre¬ 
sponding  to  an  expansion  from  to  any  pressure  P  in  equation 
(3)  is  proportional  to  the  square  root  of  the  pressure  change. 
As  a  result,  the  pressure  change  in  the  expansion  wedge  as  a 
function  of  the  angle  approximates  a  second  order  parabolic 
decay 

P  =  Pi[1_(l^)2] 

In  all  practical  cases,  the  shockwave  to  be  considered 
is  spherical.  Some  general  considerations  applicable  to  this 
case  are  in  order  here.  Because  of  the  curvature  of  the  front, 
the  angle  0,  the  pressure,  and  consequently  all  other  vari¬ 
ables  change  with  time  or  distance  from  the  charge  center;  a 
stationary  flow  can  therefore  no  longer  be  realized.  However, 
over  a  very  small  area  close  to  the  free  surface,  the  curv¬ 
ature  of  the  front  is  negligible,  and  the  stationary  flow 
pattern  illustrated  in  Figure  1  is  applicable  a_t  any  one  time 
(a  quasl-stationary  flow).  If  the  region  of  interest  extends 
an  appreciable  distance  down  from  the  surface,  the  shock  front 
and  the  expansion  wedge  become  arcs  (the  latter  only  as  an 
approximation).  Superposition  of  a  radial  flow,  parallel  to 
the  surface  at  the  surface,  as  illustrated  in  Figure  3*  makes 
the  entire  shock  front  and  the  expansion  wedge  at  the  surface 
quasi-stationary .  The  expansion  wedge  farther  down,  however, 
will  not  be  completely  quasi-stationary.  The  greater  the 
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curvature  of  the  arcs  and  the  greater  the  distance  between  the 
shock  front  and  the  expansion  wave,  the  greater  will  be  the 
deviation  from  the  quasi -stationary  flow.  For  the  geometries 
considered  in  this  paper,  the  deviation  is  not  serious  and  can 
be  neglected.  Thus,  it  becomes  possible  here  to  consider  a 
stationary  flow  pattern  even  in  the  case  of  a  spherical  shock- 
wave.  This  makes  it  feasible  to  extend  many  of  the  results 
obtained  in  the  plane  case  to  the  actual  spherical  problem  by 
an  analogy  to  the  method  of  Images  used  in  acoustics. 


Ill  THE  CRITICAL  CONDITIONS  AND  THE  BOUNDARY 
OF  THE  ANOMALOUS  REGION 

For  a  plane  shock  wave,  as  shown  in  Figure  1,  the  flow  f^ 
into  the  expansion  wedge  is  a  function  of  the  glancing  angle  0 
and  the  pressure  P.  This  is  true  because  the  propagation 
velocity  U  and  the  particle  velocity  u  depend  only  on  P.  If 
for  a  shockwave  of  given  pressure,  the  angle  0  is  decreased, 
the  absolute  magnitude  of  ^  decreases,  ai.d  eventually  f1 
may  become  sonic.  In  that  case,  the  initial  Mach  line  of  the 
expansion  wedge  must  be  perpendicular  to  the  flow  f^;  the 
flow  out  of  the  wedge  remains  supersonic.  If  the  flow  f1 
becomes  subsonic,  the  entire  Prandtl -Meyer  treatment  collapses, 
and  such  centered  expansion  at  a  corner  is  no  longer  possible. 

The  point  where  the  flow  fi  is  sonic  constitutes  the 
boundary  between  regular  and  "anomalous"  surface  reflection. 

For  the  case  of  an  infinite  plane  shockwave,  the  condition  of  a 
sonic  flow  f.^  leads  to  a  relation  between  the  shock  pressure  P 
and  the  critical  glancing  angle  ®crlt:  From  the  velocity 
triangle  in  Figure  1,  it  follows  that 
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(U 


sec©-u  cos( 


2  2  2 
r  -  u^sin^® 


(7) 


But  f^  =  c  (the  local  sound  velocity)  is  the  condition  of 
sonic  flow;  therefore, 

(U  sec«crlt  -  U  ooi»crlt)2  -  u2sin2«crlt  =  c2  (8) 

or  ®orlt  =  tan'1  {^[o2-(U-u)2]1/2j  (9) 

where  U,  u,  and  c  are  functions  of  the  shock  pressure  P. 

The  fact  that  the  first  Mach  line  of  the  expansion  is 
perpendicular  to  the  sonic  flow  f^,  leads  to  the  expression 

Xcrit  -  sin'1  if  sln8crit}  <10> 


where  is  the  angle  which  the  first  expansion  Mach  line  makes 
with  the  perpendicular  to  the  free  surface.  This  angle  becomes 
very  small  for  pressures  less  than  10,000  psi  and  can  usually 
be  neglected.  In  other  words,  the  first  Mach  line  can  be 
considered  as  being  perpendicular  to  the  free  surface. 

If  the  plane  shockwave  of  the  preceding  discussion  is 
considered  to  be  a  small  section  of  a  spherical  shockwave  at 
the  free  water  surface,  the  pressure  P  and  its  functions  U,  u, 
and  c  become  functions  of  the  radial  distance  from  the  charge 
center.  Substitution  of  the  numerical  values  given  in  the 
Appendix  transforms  equation  (9)  into 


tan  ®crit 


'1.197  »  1Q~5 
_l+5.367  *  10 


P 


P 


(11) 
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and  tan  ®crlt  =  1.820  r"*6°  -  1.353  r"1*80  Pentolite 

=  1.574  r"‘565  _  --1.695  TNT 

where  r  represents  the  radial  distance  from  the  charge  center 
expressed  in  charge  radii* . 

Thus,  for  a  spherical  wave,  there  exists  a  point  at  the 
surface  beyond  which  the  reflection  is  anomalous.  For  any 
particular  position  of  the  charge,  that  is  for  a  charge  at 
depth  d  below  the  surface,  this  point  will  be  at  a  definite 
radial  distance,  from  the  charge  center,  given  by 


?cnt  -  at1-42*  + 

=  3(1.449  + 

As  d  increases,  rcrlt 
sion  (13)  reduces  to 


.3015  rCr[l°)1//2  Pentolite 

(13) 

.4047  r^J.3)1/2  TNT 
increases,  and  for  3^12,  expres- 


rcrit  =  -224  d2,5°  Pentolite 

=  .354  32‘3°  TNT 


(14) 


This  shows  that  a  region  of  anomalous  surface  reflection 
exists  for  every  charge  depth**.  The  distance  from  the  charge 
center  to  the  beginning  of  this  region  at  the  surface, 
increases  rapidly  with  charge  depth. 


throughout  this  paper,  symbols  with  bars  over  them  will 
represent  distances  expressed  in  charge  radii. 

**Cf  course,  the  numerical  equations  apply  only  as  long  as  the 
shockwave  pressure-distance  relationships,  as  given  in  the 
Appendix,  are  valid.  According  to  reference  [C],  this  is  true 
to  relatively  large  distances. 
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Next,  consider  a  plane  shockwave  which  suddenly  interacts 
with  the  free  surface  at  a  glancing  angle  which  is  less  than 
the  critical  angle.  From  the  point  of  intersection  of  the 
shock  front  with  the  surface,  a  signal  will  originate  which 
propagates  with  the  local  sound  velocity  c  into  the  local 
particle  stream  of  velocity  u.  Since  the  glancing  angle  is 
less  than  the  critical  angle,  the  horizontal  velocity  compo¬ 
nent  of  the  signal  is  greater  than  that  of  the  shock  front. 
Thus,  the  surface  signal  overtakes  the  shock  front.  This  is 
illustrated  in  Figure  4. 

The  path  of  descent  of  the  primary  surface  signal  must 
be  the  line  along  which  the  shockwave  has  the  same  velocity 
as  the  surface  signal.  This  is  the  case  where  vector  addition 
of  c  and  u  equals  U,  or  in  other  words  where  c  equals  f^. 

This  condition  is  analogous  to  the  critical  condition  at  th<= 
free  water  surface,  so  that  the  path  of  descent  may  be  con¬ 
sidered  a  free  surface,  which  the  incident  shock  intersects 
at  the  glancing  angle  9crit*  The  shock  front  below  this  path 
is  unaffected  by  the  surface  signal,  above  this  path  it  is 
attenuated.  The  region  where  the  shock  front  is  attenuated 
constitutes  the  "anomalous"  region . 

According  to  these  considerations,  Keil  [I]  has  defined 
the  boundary  of  the  anomalous  region  as  the  curve  which 
intersects  the  radius  vector  r  from  the  charge  center  at  the 
critical  angle.  This  leads  to  a  differential  equation  for 
the  rate  of  descent  of  the  boundary: 


d?crlt 

dr 


tan  ® 


crit 


.  a-ycnt 
r  J  r 


(15) 


where  Ycrit>  expressed  in 
anomalous  region  boundary 


charge  radii,  is  the  depth  of  the 
beneath  the  free  surface.  For  a 
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FIG.  4  PROPAGATION  OF  PRIMARY  SIGNAL  ALONG  PLANE  SHOCK  FRONT 

(RESTING  POINT  OF  OBSERVATION) 
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spherical  wave,  the  critical  angle  is  a  function  of  the  radius 
vector  r,  as  shown  in  equation  (12).  Substitution  of  the 
appropriate  expressions  in  the  differential  equation  (15)  and 
integration*  lead  to  the  expressions 

Pentolite 


d-y 


crl—  =  sin  (3.033  r"*6°-.752  r_1,8°  +  .497  r*3*00 


-  k') 


TNT 


(16) 


=  sin  (2.78^  r',565-.515  f-1,695  +  .257  r-2'825  -  k') 


k*  is  a  constant  of  integration  evaluated  at  the  point  of 

critical  surface  reflection:  that  is,  at  r  =  r  y  , , 

cpit*  J crit 

must  be  0. 

In  the  preceding  derivation,  the  path  of  a  signal  origin¬ 
ating  at  the  point  of  critical  surface  reflection  was  traced. 
It  appears  reasonable  that  this  signal  is  indeed  the  first  to 
arrive  at  the  shock  front  at  a  given  depth  yj  however,  no 
proof  is  yet  at  hand  that  a  signal  originating  at  the  surface 
prior  to  the  point  of  critical  reflection  must  arrive  at  the 
shock  front  at  a  later  time. 


IV  THE  ANOMALOUS  REGION 

The  discussions  of  the  previous  section  lead  to  a  defin¬ 
ite  picture  of  the  shockwave  in  the  anomalous  region:  a  spher¬ 
ical  shock  front  Invades  the  undisturbed  water.  An  expansion 
wave  originating  from  the  free  surface  has  overtaken  the 
incident  shock  front  near  the  surface,  and  its  head  intersects 
the  shock  front  at  the  known  depth  ycrj_t«  The  shock  is  atten¬ 
uated  when  overtaken  by  this  rarefaction  wave.  The  shock 

*W.  Hostow,  in  a  private  communication,  pointed  out  that  the 
differential  equation  (15) >  due  to  Keil  [I],  can  be  integrated 
analytically. 
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pressure  will  decrease  from  the  free  water  value  at  ycrit  to 
a  value  P  at  the  surface,  in  a  manner  depending  on  the 

o 

pressure  profile  of  the  rarefaction  wave. 

At  this  point  the  assumption  is  made  that  the  conditions 
in  and  immediately  below  the  anomalous  region  can  be  repre¬ 
sented  approximately  by  a  quasi-stationary  flow  pattern 
together  with  a  nonstationary  rarefaction  wave.  Figure  5 
illustrates  such  a  flow  pattern  for  a  plane  shockwave.  A  flow, 
-Up  sec  0,  which  makes  the  undistorted  shock  front  stationary, 
is  superposed.  The  head  of  the  nonstationary  rarefaction  wave 
intersects  the  shock  front  at  the  boundary  of  the  anomalous 
region,  that  is,  at  the  depth  ycr^t*  The  point  of  intersec¬ 
tion  of  the  actual  shock  front  with  the  surface  must  coincide 
with  the  stationary  rear  of  the  nonstationary  rarefaction 
wave.  The  remaining  pressure  must  be  released  at  this  point; 
this  is  accomplished  by  a  stationary  centered  expansion  wave, 
the  Prandtl-Meyer  corner  flow. 

The  pressure  P  at  the  rear  of  the  nonstationary  rare- 

o 

faction  wave  can  be  found  from  the  condition  that  the  flow 
should  be  stationary  there.  This  leads  to  the  following 
procedure:  The  flow  into  the  nonstationary  wave  is  split 
into  components  parallel  and  perpendicular  tc  the  wave. 


Parallel  component  = 

(U.-,  sec©  -  u„  cos®)  sinY  -  u„  sin®  cosy 


(17) 


Perpendicular  component  = 


(Up  sec®  -  Up  cos®)  cosy 


(18) 


where  y  ■  (®crit  ”  '  The  wave  accelerates  the  perpendicular 
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component  of  the  flow.  For  a  plane  wave,  the  acceleration 
depends  only  on  the  overall  pressure  change.  At  the  rear  of 
the  wave,  the  perpendicular  component  of  the  flow  must  be 
equal  to  the  local  sound  velocity.  Thus, 

rPF  v 

cs  =  (Up  sec®  -  Up  cos®)  cosy  +  \  —  dP  (19) 

Ps 

The  last  term,  the  Riemann  function,  represents  the  velocity 
increment  which  the  nonstationary  wave  imparts  to  the  flow. 

The  expression  can  be  solved  for  the  pressure  P  at  the 
beginning  of  the  stationary  expansion  centered  at  the  surface. 

Actually,  the  shock  front  in  the  anomalous  region  is  bent 
back  and,  because  of  its  lower  peak  pressure  there,  propagates 
at  a  lower  velocity  than  the  unattenuated  shock.  As  a  rough 
approximation,  it  can  be  assumed  that  these  two  effects  cancel 
as  far  as  the  horizontal  propagation  velocity  of  the  wave  is 
concerned .  In  other  words , 


Up  s^c®  =  Uy  sec0y  =  Ug  sec0s 


(20) 


where  pf  represents  the  glancing  angle  of  the  attenuated  front, 
and  the  subscripts  s  and  y  refer  to  the  surface  and  to  any 
depth  y  in  the  anomalous  region  respectively  (see  Figure  5). 
Equation  (20)  states  that  the  superposed  flow  Up  sec®  makes 
the  entire  shock  front  stationary  with  respect  to  the  surface. 
Substitution  of  the  numerical  values  from  the  Appendix  into 
this  equation  leads  directly  to  a  relation  between  0  and  P  , 

v  w 

the  peak  shock  pressure  at  any  depth  y: 


P  =  1.863  *  105 

V 


Sr 


(1  +  5.367  *  10'  Pp)  cosft 


cos® 


*  -  1 


(21) 
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It  can  be  shown  that  rotation  is  introduced  into  a  stream 
on  passing  through  a  bent  shockwave  [0].  In  the  present  case, 
this  effect  causes  the  rarefaction  wave  to  bend  back  slightly 
close  to  the  surface.  At  the  relatively  low  pressure  of 
interest  here,  this  bending  is  small  and  can  be  neglected. 

Next,  equation  ( 19 )  is  evaluated  for  the  quasi-station- 
ary  case  of  a  plane  wave,  whose  parameters  change  with  dis¬ 
tance  from  the  origin  like  those  of  a  spherical  shockwave . 

This  leads  to  the  expression 


't**- 


1.570  *  105  dr"1,Go  -  3.63  *  103  r‘1,20  Pentolite 

(22) 

1.360  *  105  dr-1  *5^5  _  3.3Q,  io3  f"1,13  TNT 


For  the  actual  spherical  case,  an  analogy  to  the  method 
of  images  permits  arcs  to  be  substituted  for  the  plane  wave 
configuration  in  Figure  5.  Account  must  be  taken  of  the 
spherical  divergence  term  arising  from  the  equation  of  conti¬ 
nuity.  This  term,  which  must  appear  on  the  right  hand  side  of 
equation  (19) ,  is  approximated  by 


.'“total 

P  dx 
o 


(r  »  % 


*  total ^ 


(23) 


where  x  is  the  distance  from  the  back  of  the  nonstationary 
wave  perpendicular  to  its  direction  of  propagation,  xtQtal  is 
its  total  width  and  pQ  is  the  density  of  the  undisturbed  water. 
(23)  is  conveniently  expressed  as  a  correction  term 


&P 


q  x 


total  F 


spherical 


(24) 


which  must  be  added  to  the  incident  surface  pressure  of  the 
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plane  case,  equation  (22).  The  symbol  q  in  equation  (24) 

represents  a  fraction,  which  depends  on  the  pressure  profile 

of  the  rarefaction  wave.  Its  value,  based  on  experimental 

observations,  appears  to  be  approximately  .85. 

The  pressure  profile  of  the  rarefaction  wave  is  related  by 

a  geometric  argument  to  the  pressure-depth  variation  along  the 

shock  front  in  the  anomalous  region.  The  latter  is  related  to 

the  angle  0  of  the  shock  front  in  the  anomalous  region  by  the 
«/ 

stationary  flow  assumption  and  the  Rankine-Hugoniot  conditions 
[see  equation  (21)].  From  various  qualitative  arguments,  the 
pressure,  and  therefore  also  the  angle  0  ,  Is  expected  to 
change  most  rapidly  with  depth  near  the  surface  and  most  slowly 
at  the  lower  edge  of  the  anomalous  region.  This  suggests  a 
power  law  for  the  angle  0  as  given  by  the  expression 

«7 

0y  =  (0a  -  ®) (1  -  y)n  +  ®  (25) 

where  f  =  y/7crit  Is  the  "relative  depth"  in  the  anomalous 
region,  and  n  is  a  constant,  whose  value  from  actual  experi¬ 
ments  appears  to  be  approximately  4. 

Knowledge  of  the  angle  0  as  a  function  of  depth  permits 

y 

a  rough  estimate  of  the  total  lag  of  the  actual  shock  front  at 
the  surface  behind  one  which  is  not  attenuated  by  surface  sig¬ 
nals.  This  is 
wave  case: 

^total 

which  becomes 

Ltotal 

when  $  =  (0S  - 
behind  that  at 


given  by  the  following  equation  for  the  plane 


■J! 


crit 


tan  0  dy  -  yCI>u 


tan0 


yorlttan  5 

(n  +  1) 


(26) 


(27) 


0)  is  relatively  small.  The  lag  at  the  surface 
any  depth  y  in  the  anomalous  region  is  similarly 
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given  by 


T'crlt  t3nS 

(n  +  1  j 


[1 


U-y)n+1l 


(28) 


These  equations,  which  were  derived  for  plane  waves,  are  also 
approximately  correct  for  the  spherical  case. 

The  pressure  of  the  nonstationary  rarefaction  wave  is 
constant  alone  a  line  for  which  x,  defined  in  connection  with 
equation  (23),  is  constant.  Thus,  the  pressure  profile  of 
the  rarefaction  wave  follows  from  the  profile  of  the  shock 
front  in  the  anomalous  region  [given  by  equations  (21)  and 
(25)  as  discussed  previously]  if  the  corresponding  values 
of  x  and  y  along  the  shock  front  are  known.  This  relation 
between  x  and  y  follows  from  geometric  considerations  and  is 
given  by  the  approximate  expression 


x  « 


tan#T  [l-(l-y)n+1]  +  y  +  sirvy] 


n  + 


-  y 


2  ycrit 


2  r 


tan§ 
n  +  1 


*  +  ycrlt  .  ycrit 

- - -  +  sinY  -  -g-y- 


(29) 


where  x  =  x/x^0tai  measures  the  relative  distance  through  the 
nonstationary  wave.  Note  that  x  ■  0  at  the  stationary  end  of 
the  wave,  and  3c  =  1  at  the  head  of  the  wave. 

The  flow  through  a  stationary  expansion  wave  was  solved 
for  the  plane  case  in  Section  II  of  this  report.  Thus,  in 
equation  (3),  Ps  must  be  substituted  for  P^,  and  the  total 
velocity  of  the  flow  into  the  wedge  is  given  by 


(1  +  8.449  *  10~6Pj2  + 

b 


(1+2.463 x  io"6P„)  sinYl2 


cos® 


(20) 


This  provides  everything  to  calculate  the  total  deflection 
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of  the  Mach  lines  at  the  Prandtl-Meyer  corner.  From  this  and 
from  the  definition  of  7,  stated  In  conjunction  with  equations 
(17)  and  (18),  It  follows  that  y  ,  the  angle  which  the  rear  of 
the  Prandtl-Meyer  expansion  wedge  makes  with  the  perpendicular 
with  the  surface.  Is  given  by 

Y  -  Y  +  (31) 

It  appears  worthwhile  here  to  make  a  few  remarks  on  the 
assumption  of  stationary  flow,  which  was  Introduced  In  this 
section.  Actually,  the  plane  wave  case  consists  of  a  plane 
incident  shock  front  and  an  approximately  cylindrical  expan¬ 
sion  wave,  which  is  attenuating  the  incident  shock.  The 
attenuated  part  of  the  shock  front  is  traveling  at  a  slower 
rate  than  the  unaffected  part  and  thus  eventually  bends  back 
near  the  surface.  Strictly  speaking,  it  is  only  possible  to 
make  the  unattenuated  part  of  the  shock  front  stationary  at 
any  one  time;  however,  the  theory  derived  In  this  section 
assumes  a  stationary  flow  through  the  entire  shockwave  config¬ 
uration.  The  approximate  character  of  this  treatment  shows 
up  In  nonstationary  effects  at  the  surface:  the  superposed 
flow  through  the  nonstationary  rarefaction  wave  Is  not  quite 
perpendicular  to  the  wave  front.  The  flow  through  the 
stationary  rear  of  this  wave  Is  thus  slightly  greater  than 
sonic,  a  condition  which  tends  to  Increase  the  Incident 
pressure  at  the  surface  point. 
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V  THE  CONSTRUCTION  OF  PRESSURE -TIME  HISTORIES 

If  the  theory  developed  in  the  preceding  sections  is  to 
be  of  practical  significance,  it  must  predict  satisfactorily 
the  pressure-time  history  of  a  shockwave  at  a  distance  r  from 
the  charge  center  and  at  a  depth  y  below  the  free  surface. 

A  procedure  which  yields  approximate  pressure-time  histories 
is  outlined  in  this  section. 

To  begin  with,  the  quantities  r,  d,  6,  ®crlt>  ycrlt> 
and  P0  for  the  plane  case  are  calculated  from  the  known  charge 
weight  and  geometry  of  the  problem  by  means  of  the  appropriate 
equations  In  the  text.  A  scaled  diagram  of  the  shock  front 
is  then  drawn,  the  bending  back  effect  near  the  surface  being 
neglected  at  first.  At  the  critical  depth,  ycrlt,  on  the 
shock  front,  a  line  is  drawn  which  makes  the  angle  with 

the  radius  vector  F  from  the  charge  center.  The  line  is 
extended  back  until  It  Intersects  above  the  free  surface  a 
line  passing  through  the  charge  center  and  perpendicular  to 
the  surface.  The  point  of  intersection,  A  in  Figure  6,  which 
in  a  sense  represents  an  "image"  point  for  the  nonstationary 
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rarefaction  wave,  serves  as  a  center  for  two  concentric  cir¬ 
cles  passing  through  the  critical  point  on  the  shock  front 
and  through  the  surface  point  of  the  shock  front  respectively. 
The  distance  between  these  circles  represents  a  first  approx¬ 
imation  to  the  width,  x.  .  ,  of  the  nonstationary  wave.  It 

is  now  possible  to  calculate  APSpherical  and  Ps  sP^er“ 

ical  cane,  and  hence  also  $s,  and  LtotQl.  If  the 

latter  quantity  is  appreciable,  it  must  be  included  in  x.  , 

to  uSl 

and  values  for  P  ,  QS  ,  etc.,  may  have  to  be  recalculated. 

s  s 

Finally,  y  is  calculated,  and  the  corresponding  "image'1  point, 
B  in  Figure  6,  is  found  by  drawing  a  line  through  the  correc¬ 
ted  surface  point  of  the  incident  shock,  which  makes  the  angle 
Y  with  the  free  surface.  The  point  of  intersection  of  this 
line  with  the  perpendicular  to  the  surface  above  the  charge 
serves  as  the  center  of  a  circular  arc  representing  the  end  of 
the  centered  expansion.  The  final  diagram  is  as  indicated  in 
Figure  6. 

Next,  a  path  through  the  shockwave  configuration  is  con¬ 
sidered  which  corresponds  approximately  to  the  way  the  shock- 
wave  passes  a  gauge.  All  distances  are  then  converted  to  a 
time  scale.  It  is  sufficient  to  use  the  sound  velocity  in 
undisturbed  water  for  this  purpose.  The  peak  shock  pressure 
and  the  pressure  decay  due  to  the  nonstationary  and  the  sta¬ 
tionary  centered  wave  are  calculated  from  the  appropriate 
equations  in  the  text.  The  results  are  superposed  on  the 
pressure-time  decay  curve  of  the  free  water  shockwave  as  given 
in  equation  (a4)  of  the  Appendix,  and  the  resulting  curves  are 
smoothed  out.  The  approximate  pressure-time  histories  thus 
drawn  represent  the  results  obtainable  with  the  present  theory. 

In  conclusion,  it  may  be  worthwhile  to  discuss  the  range 
of  applicability  of  the  present  work.  In  general,  it  is 
possible  to  calculate  pressure-time  histories  for  charge 
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depths  of  3  charge  radii  or  more,  at  distances  where  the  free 
water  shockwave  peak  pressure  is  10,000  psi  or  less.  In  this 
case,  the  equations  are  all  in  analytical  form,  and  agreement 
with  experiment*  appears  satisfactory.  It  should  be  noted 
that  the  present  theory  was  derived  for  the  ideal  case  of  a 
semi-infinite  homogeneous  medium  and  a  perfectly  plane  free 
boundary.  However,  under  actual  experimental  conditions, 
other  factors  enter  -  such  as  surface  roughness,  temperature 
and  density  gradients  in  the  undisturbed  medium  [G],  the 
presence  of  a  solid  bottom,  etc.  Their  effects  depend  on  the 
size  and  locale  of  the  particular  explosion  and  must  be 
investigated  separately.  Finally,  it  should  be  pointed  out 
that  the  present  results  do  not  apply  to  oharge  depths  of  less 
than  2  charge  radii.  Here  the  assumptions  of  the  theory  are 
no  longer  satisfactory,  and  in  addition,  the  very  difficult 
problem  of  an  air-backed  explosion  arises.  It  is  felt  that 
for  these  very  shallow  charge  depths,  more  work,  especially 
experimental,  is  necessary. 


*  See  References  [E]and  RJ  . 
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APPENDIX 

PREE  WATER  SHOCKWAVE  PARAMETERS 

For  a  shock  front  moving  through  a  fluid  medium,  the 
hydrodynamic  and  thermodynamic  transitions  are  specified  by 
the  Rankine-Hugoniot  relations;  for  an  expansion  wave,  on  the 
other  hand,  isentropic  change-,  are  applicable.  If  the  initial 
state  of  the  medium  is  known,  and  if  sufficient  P-V-T  data  is 
available,  both  cases  can  be  expressed  completely  in  terms  of 
one  independent  variable  -  for  example  the  pressure  change. 

For  the  relatively  low  pressure  changes  applicable  to  this 
report,  the  isentropic  and  the  Hugoniot  curves  for  water  are 
almost  identical  and  differ  but  little  from  the  isothermal. 

The  following  Rankine-Hugoniot  parameters  for  fresh  water, 
initially  at  rest  and  at  20°C  and  1  atm.,  have  been  calculated 
with  the  aid  of  references  [A,H,J]  and  are  quite  satisfactory 
up  to  pressures  of  30,000  psl .  For  much  higher  pressures,  the 
numerical  values  given  in  reference  [P]  should  be  used.  The 
following  relations  hold  for  this  low  pressure  region: 

U  -  cQ(l  +  5-367  * 10"6  P) 
u  =  cQ(2.89  *  10-6  P) 

U  -  u  =  c  (1  +  2.463  x  IQ-6  P)  (A1) 

0 

c  =  c0(l  +  8.449  *  10"°  P) 

where  P  is  the  excess  pressure  in  psl,  U  is  the  propagation 
velocity  of  the  shockwave,  and  u  and  c  are  the  particle 
velocity  and  local  sound  velocity  behind  the  shock  front. 
cQ,  the  sound  velocity  in  the  undisturbed  water  (20°C  and 
1  atm.)>  is  1483  meters/second. 

An  equation  of  state  for  water,  applicable  to  the  true 
as  well  as  the  Hugoniot  adiabatic  up  to  a  pressure  of  30,000 
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psi,  is  given  by  the  expression 

V  =  1.00177  (1-3.04  *  10"6P  +  2.3  *  10"11?2)  (A2) 


where  V  is  the  specific  volume  of  the  water,  and  an  initial 
state  of  20°C  and  1  atm.  is  assumed. 

The  experimental  free  water  shockwave  peak  pressure  data 
are  conveniently  expressed  in  terms  of  the  initial  charge 
radius  For  pressures  below  30,000  psi 


Pp  =  2.77 *  105  r-^*^  Pentolite,  p  =  I.63  gm/cc  [F] 
Pp  =  2.07  *  105  r_1  -13  TNT,  p  =  1.55  gm/cc  [D] 


(A3) 


where  Pp  is  the  free  water  peak  pressure  in  psi,  r  is  the 
radial  distance  from  the  center  of  the  charge  expressed  in 
charge  radii,  and  p  is  the  loading  density. 

The  free  water  pressure  decay  behind  the  shock  front, 
observed  at  a  fixed  point,  is  given  by  the  expression 

P  =  PFe"t//e  Pp  ^  15,000  psi 

where 

0  =  31.1  WV3  Pentolite,  /O  =  I.63  gm/cc  [F](a4) 

0  =  41.9  W1/3  r’18  TNT,  p  =  1.55  gm/cc  [D] 

Here  P  is  in  the  same  units  as  Pp,  t  is  the  time  after  shock 
front  arrival  in  microseconds,  0  is  the  time  constant  in 
microseconds,  and  W  is  the  charge  weight  in  lbs.  This  formula 
is  valid  only  up  to  t  «  0. 
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